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Abstract 

An upper limit to the critical density for the transition to the deconfined phase, at zero temper- 
ature, has been evaluated by analyzing the behavior of the gluon condensate in nuclear matter. Due 
to the non linear baryon density effects, the upper limit to the critical density,pc turns out about nine 
times the saturation density, po for the value of the gluon condensate in vacuum < (Qs/7r)G^ >= 0.012 
GeV'^. For neutron matter pc — 8.5po- The dependence of the critical density on the value of the 
gluon condensate in vacuum is studied. 



1 Introduction 



The future experimental results both in relativistic heavy ions collisions and in astrophysics will give a 
definite answer on the observation of new states of matter under extreme conditions. Our understanding 
of particles and nuclear phenomena suggests that at finite temperature there is a phase transition to 
deconfined gluons and massless quarks. QCD lattice calculations H strongly support the conclusion that 
at a critical temperature, Tc, a transition to the quark-gluon plasma phase occurs and Tc is equal to T^^ 
the critical temperature for the restoration of the chiral symmetry. More recently it has been proposed 
that, at high density and low temperature, QCD exhibits a new phase due to a diquark condensation, 
called color superconductivity, with interesting astrophysical effects [Q- The understanding of the whole 
phase diagram of the theory requires an analysis of the dependence of the order parameters on the 
temperature and on the density. Unfortunately lattice calculations at large density are problematical 
although recently two methods have been proposed to overcome, at least for not too large values of the 
chemical potential, the technical difficulties that arise in lattice simulations 3 . 

The regime, of low T and high density, is particularly interesting not only from the theoretical point 
of view but also because it has strong implications in astrophysical systems such as the neutron stars. 
For example, the mass and the relation between the mass and the radius of a neutron star depends on 
the state of quark matter in the core ( if any) and on the transition from nuclear matter to quark matter 
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The value of the critical density at which this possible transition occurs inside the neutron star has 
been determined, up to now, by comparing and matching the two Equations of State (EoS), one for 
the nucleon matter and one for the deconfined quark matter, which are derived within two different 
frameworks. In this procedure a first order phase transition is assumed, and the Gibbs construction, or 
a modification of it, is employed. 

For the nucleon EoS, either microscopic or phenomenological approaches have been followed. In the 
microscopic approach, the nucleonic EoS is derived from a realistic nucleon-nucleon (NN) interaction, 
which fits the NN phase shifts and deuteron data, by solving accurately the many-body problem for 
infinite nuclear matter. The main requirement is the reproduction of the empirical saturation point of 
symmetric nuclear matter, as extracted from the Weizacker-like mass formula for finite nuclei, and the 
compressibility at saturation as extracted from the analysis of monopole vibrations in nuclei. As a rule, 
three-body forces have to supplement the two-body forces in order to fulfill these requirements. In the 
purely phenomenological approach, the same requirements are used to fix or to constrain the parameters 
of the model 

For the quark phase, only simple models can be used, since no "a6 mitio" QCD calculations at high 
baryon density is available. Both MIT bag model and the NJL model have been used. Due to the 
uncertainties inherent to the models, the precise location of the phase transition is not well known. 

Moreover another constraint on the value of the critical density, is the phenomenological observation 
that in heavy ion collisions at intermediate energy ( lOMeV/A < E/A < 200MeV/A ) no evidence of a 
transition to the quark-gluon plasma phase has been found. Indeed, all microscopic simulations that are 
able to reproduce a great variety of experimental data, do not need the introduction of such a transition. 
In these simulations the calculated nucleon density can reach values which are at least 2-3 times larger 
than the saturation density pQ. One can, therefore, conclude that symmetric or nearly symmetric nuclear 
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matter at few MeV of temperature does not exhibit any phase transition to the deconfined matter up 
to this baryon density. It has to be noticed that the phase transition in symmetric matter can occur at 
a different baryon density than in neutron star matter, where nuclear matter is closer to neutron matter 
than to symmetric nuclear matter. 

This constraint coming from heavy-ion physics appears as an independent one, that should be fulfilled 
by any theory or model of deconfinement. Indeed, the quark matter model mentioned above seem to 
have, in some cases, serious difficulties to fulfill the constraint (the transition occurs at too low density, 
or it does not occur at all), although it produces " reasonable " results for neutron stars, where the 
transition does occur. According to the previous discussion, it would be, therefore, desirable to formulate 
an independent and differently founded estimate of the transition point (or the possible mixed region). 

In this Letter we shall address this problem by discussing the dependence of the QCD gluon condensate 

< — (0)G'^'"^(0) >= ^ (1) 

TT 

on the density p (at zero temperature), as a signal of confinement-deconfinement transition. In fact the 
trace anomaly decreases with temperature but does not vanish above the critical temperature and 
then the gluon condensate is not, by itself, an order parameter for the transition. However, the vanishing 
of the gluon condensate at large density gives an upper limit to the critical density for the following 
reasons: 

1) According to Dosch |Hj the string tension a is given by 

a^M' J drdrf C^'j"^' ) (2) 

where a is the correlation length for the gluon field defined by 

< ^ G;i,(a;)G^-(0) >= M^fix^/a") (3) 

TT 

By Equation |(2Jl, it turns out that the string tension in proportional to 

a ~ M^a^ (4) 

This implies that the behaviour of the string tension at finite density depends not only on the depletion 
of the gluon condensate, M**, but also on the decreasing of the gluon correlation length at large density. 
Therefore the vanishing of the gluon condensate can give only an upper limit to the critical density. 

2) The gluon condensate is the sum of two independent terms, < > and < >, see for instance 
the review in jH). At finite temperature, < > vanishes at the deconfinement transition, while < > 
remains finite even above the critical point and one can expect an analogous behavior at finite density 
[TUl im [T^ . Again, this is an indication that the full gluon condensate provides an upper limit to the 
critical density. 

3) As discussed by Hatta and Fukushima JSIj the vanishing of the glueball can describe the critical 
phenomena at the deconfinement transition. But the glueball mass is usually considered proportional to 
the trace anomaly, 

Mg ^< — (0)G^""(0) >i/^= M (5) 

TT 
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which persists at any temperature. Then there is an apparent contraddiction solved by considering that 
only a fraction of the glueball mass comes from the trace anomaly ( see ^3] ) and that, as previously 
commented, only the "electric" glueballs (those which contains time-like links) become massless at the 
phase transition while the "magnetic" glueballs remain massive. Then the vanishing of the glueball mass 
occurs at a smaller density than the estimate based on the gluon condensate. 

According to the above issues, we shall evaluate the upper limit to the critical density by analyzing 
the dependence of the gluon condensate on the baryon density. In Section 2 we shall give an heuristic 
argument for the occurence of the transition, in Section 3 the known results about the behavior of the 
gluon condensate in nuclear matter are recalled and we discuss a phenomenological consequence of this 
behavior; Section 4 is devoted to the calculation of the critical density including the non linear effects 
and Section 5 to the conclusions. 



2 TOY CALCULATION 

The glueball mass is related to the gluon condensate. Then the first heuristic argument about the critical 
density can be obtained by looking at the behavior of the glueball mass as a function of the baryon 
density. In this Section we shall consider a simple calculation to give some very preliminary indications 
of the deconfinement transition at large density. 

Let us consider the color potential given by a linear rising term plus the perturbative contribution 

V{r)=ar--- (6) 
3 r 

where a is the effective string tension, a is an effective constant (independent on the scale) and the factor 
4/3 is due to the color degrees of freedom. Let us assume that the glueball of total energy i? is a bound 
state of two massless gluons: 

4 a 

E^2p + ar~-- (7) 
3 r 

From the uncertainty principle pr > 1 and then, by minimizing E with respect to r, it turns out that 
otq = (2 — 4a/3) and E{r — vq) = Eq ~ 2arQ. For (almost) reaUstic values a ~ 0.3 GeV"^ and a ~ 0.1, 
i?o — 1.5 GeV not far from the present glueball mass limit. 

This simple picture is modified by baryon density effects ( which will be discussed later) and in the 
medium the potential in Eq. jnj is replaced by the screened potential ^1] 

1 - c"'"" 4 a 
/xr 3 r 

where n is the screening parameter and the screened potential has the asymptotic limit Vs{oo) = cr/zi. By 
repeating the previous calculations, it turns out that rp = rf){p) , and the value of the minimum energy 
is 

rr r / 1 \ 1 

(9) 

which shows that there is a critical value of the screening parameter at which the bound state disap- 
pears. 
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To evaluate the dependence of Eq on the baryon density p (and then the critical density) we need to 
know its relation with the screening parameter. By assuming a Debye screening due to a 2 flavor massless 
quarks gas of density pq = 3/5 ^Hl j then 

, 12a / \ ^ 

By means of Eqs. 10) and it turns out that the critical value where the bound state disappears is 
Pc — 0.24 GeV that corresponds to a baryon density p ~ 3po where po ~ 0.16 fm~^ is the saturation 
density. 

The previous argument represents only a qualitative way to show that the glueball disappears from the 
spectrum at some critical density and does not pretend to make any quantitative statement. For example, 
if one fixes the value of the parameter a in the potential V{r) to fit the glueball mass nig ~ 1.7 GeV, 
then it turns out that the critical density is about five times the saturation density. 



3 GLUON CONDENSATE IN NUCLEAR MATTER 

In this Section we recall the main results concerning the behavior of the gluon condensate in nuclear 
matter (see ^7]). In vacuum, at T = 0, the estimates of the gluon condensate ^Hl indicate the range 

0.003 GeV^ < (< ^G°^(0)G^''°(0) >) < 0.021 GeV^ (11) 

By defining the gluon condensate at zero temperature and at finite density in nuclear matter 

gip)^ < Af |^G;^(0)G^-(0)|A/ >, (12) 

TT 

the non linear density effects are taken into account by considering |17| 

gip) = .g(0) - Y + e(p) - ^ m, < N\q^q,\N >) (13) 

where m is the nuclcon mass, rrii are the current quark masses, < N\q^qi\N > is the fermionic condensate 
in the nuclcon in the medium, the gluon condensate in vacuum is indicated with g(0), e{p) is the binding 
energy per nucleon in the medium and the sum of the fermionic condensates is extended only to the three 
light flavors [T71IB] . 

According to previous equation the gluon condensate at finite density decreases by increasing p and, 
despite this result refers to nuclear matter, one can ask if some phenomenological indication of this 
depletion comes out, for instance, by comparing light and heavy nuclei. To answer to this question one 
can consider the behavior of the deep inelastic scattering structure function, F^, per nucleon in a nucleus. 

For very small values of the scaling variables x the experimental data on the ratio of the structure 
functions for two different nuclei, F^/F.^, shows a depletion (shadowing). 

In the small x region the structure function is described by the pomeron exchange |2()j and, indeed, 
for a; — > the structure functions are related to the quark-pomeron effective coupling in the target (T), 
(3t, 

F^{x)-^Bt (14) 
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and the ratio of the structure functions per nucleon in the nucleus A with respect to nucleus B is PI] 

i^2(x)^'~» /3b ^ 

On the other hand, in the model proposed by Landshoff and Naehtmann the pomeron exchange is 
described by nonperturbative gluons and the quark pomeron coupling is related to the gluon condensate 
according to relation [521 pi] 



(3T^a^<T\ — Gl^iO)G^"'''iO)\T> (16) 

TT 



In this analysis one can consider the correlation length a as independent of the target, because it is 
related to the pomeron coupling to off-shell quarks |23| |U therefore one obtains 

F2{x)^ Pa _ <^|^G°,(0)G'^-(0)|A> 

F2{x)B^^'-o Pb <i?|^G-,(0)GA'-(0)|i3> ^ ' 

Then the last ratio can be written as 

where pA and ps are the nuclear densities of A and B. For example the experimental data on the 
structure function at small x give a ratio F2{x)^'^ / F2{x)^ of about 0.9 at a; = 0.007 ^5] . 

One can conclude that, according to the model in [22], there is at least an indication that at zero 
temperature the gluon condensate decreases in going from light to heavy nuclei. Of course the calculation 
of the numerical value of the shadowing suppression is beyond the scope of this paper. Indeed we are 
considering the nuclear matter case and not finite nuclei and moreover a quantitative prediction of the low 
X behavior of the structure functions, even in the case of the free nucleon case, needs a deep understanding 
of the transition between the perturbative and the nonperturbative regimes in QCD. 



4 NON LINEAR EFFECTS 

The indications provided in the previous Section suggest that Eq. (|13|l could be a good starting point to 
evaluate the critical density by introducing the non linear effects in e {p). This requires the analysis of 
the dynamical description of nuclear matter. 

At low baryon density and zero temperature nuclear matter can be described as a dilute gas of 
nucleons. The residual nucleon- nucleon interaction is mainly mediated by mesons. In this regime it is a 
reliable approximation to model nuclear matter as a gas of nucleons interacting through a static potential. 
The latter can be extracted from the phenomenological analysis of nucleon-nucleon scattering data. Once 
the static nucleon-nucelon interaction is given, the many-body problem can be accurately solved and the 
energy density as a function of baryon density can be calculated ( see eg. j26|). 

Only at very low density, much lower than the saturation density, the energy density follows a linear 
behavior, since different correlations start to develop in the system as the density increases. Therefore, 
the mean field approximation turns out to be inadequate. In particular, Pauli principle strongly modifies 
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the two-nucleon scattering process in the medium and this effect is density dependent. Other many- 
body effects, lifce the momentum dependence of the single particle potential and three-body correlations 
contribute to the nonlinear terms. Much work has been done along these lines, and the EoS of nuclear 
matter, up at least to saturation density, can be considered well established tZl)| . 

Near saturation density a non- negligible contribution of three-body forces and/or relativistic effects 
appear necessary to get an accurate saturation point. Their contributions is anyhow at few percents 
levels and do not affect the main trend of the energy density and of the corresponding nonlinear terms. 

At higher density the calculations contain a certain degree of extrapolation, since both two-nucleon 
and three-nucleon forces must be extended beyond the values of the relative momenta where they have 
been phenomenologically checked. However, up to density 2-3 times larger than the saturation density 
the nuclear matter EoS obtained along these lines can be considered still reliable. 

The calculation of e (p) following the previous indications has been performed in ^7\. The results 
are displayed in Table 1 of [27| for two different potentials which include three-body forces both for 
symmetric and neutron matter and we do not report them here. From these values and from Eq. H13|l . 
one obtains an approximated dependence of the gluon condensate on the baryon density. The quark 
condensate contribution, which as already noticed is restricted to the light flavors, is small because of the 
size of the quark masses with respect to the nucleon mass which also appears in Eq. . Our calculation 
is performed in the chiral limit with to„ = = and we have included the strange contribution taking 
TTis = 0.1 GeV and < A^|ss|7V >= 1, according to the estimates quoted in ^2] and the overall effect of 
the strange condensate is only few percent on the estimate of g{p)- In Eq. H13|l a possible non- linear 
contribution due to the strange quark condensate has been neglected because the included linear term 
is small and one expects that this further correction should be even smaller since the strange meson 
exchange plays a negligible role in N-N interaction, as discussed in TT . 

According to the magnitude of the condensate in vacuum ( see Eq. I^lj) one finds different values 
of the upper limit to the critical density. The behavior of the condensate as a function of the density is 
depicted in Figure 1 where g{p) is plotted for three representative values of g(0) and with e {p) determined 
in in the case of symmetric matter for the two potentials AV14 and Paris. Figure 1 remarkably shows 
the importance of the nonlinear effects. Indeed the curves that would be obtained by neglecting e{p), 
correspond approximately to the straight lines tangent to the curves in Fig. 1 at the lowest value of p. 
These curves would yield a much larger critical density (above 2.5 fm^^). The nonlinear effects due to 
e{p) reduce the critical density of almost a factor two. In the case of neutron matter the behavior is 
analogous to the one shown in Fig. 1 and, as expected, a smaller critical density is found. 

In Figure 2 the critical density is plotted as a function of the magnitude of the gluon condensate in 
vacuum both for symmetric and neutron matter. In each case we have displayed two curves corresponding 
again to the two potentials considered in j27| . Our upper limits are compatible with the values of the 
critical density suggested in 



5 CONCLUSIONS 

In summary, we have considered a way of estimating the upper limit to the critical density for the 
transition to the deconfincd phase at zero temperature by taking into account, in the determination of 
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the gluon condensate in nuclear matter, the effects of the binding energy per nucleon which have already 
been determined by means of the EoS for nuclear matter. Points 1) - 3) in the Introduction clearly 
explain that the gluon condensate is not an order parameter and that the present approach can only 
give an upper limit to the critical density. On the other hand, our analysis shows that there are good 
indications that Eq. H13|l . once the nonlinear effects have been included, can reasonably describe the 
gluon condensate dependence on the baryon density at zero temperature. 

The quantitative analysis at large baryon density gives for the (central) value of the gluon condensate 
in vacuum g{0) = 0.012 GeV^, an upper limit pc — Being an upper limit, this value is compatible 

with the typical values obtained by using phcnomenological models for the quark matter and the Gibbs 
construction OE]. A lower value of the vacuum gluon condensate, g{0) w 0.008 — 0.01, which is well 
within the theoretical uncertainty, would produce a critical density closer to these phcnomenological 
estimates. Despite these uncertainties, the result appears highly non-trivial, since the method employed 
here is totally different. The decreasing of the critical density for asymmetric matter ( pc — 8.5po)i is less 
than expected on the basis of detailed dynamical calculations but it is remarkable that this effect is 
nevertheless obtained in a totally different framework. 
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Figure 1: g{p) in GeV^ vs. the density p in fm"^ for .g(0) = .003 GeV^ (a), g{0) = .012 GeV^ (b), 
g{Q) — .021 GeV^ (c). Dashed and sohd hnes correspond to the two different determinations of e{p) in 
|27j . obtained respectively with the AV14 potential and the Paris potential, for symmetric matter (see 
text). 
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Figure 2: The critical density pc in fm^'' as a function of the condensate in vacuum g{0) in GeV^. 
Curves a.l and a. 2 correspond to the two determinations of e{p), obtained respectively with the AV14 
potential and the Paris potential, for symmetric matter as in Figure 1, and b.l and b.2 the same for 
neutron matter. 
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